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Abstract 

The radii of nucleon distribution, bulk density, and neutron skin in nuclei beyond the /3-stability 
line are studied within the direct variational method. We evaluate the partial equation of state of 
finite nuclei and demonstrate that the bulk density decreases beyond the beta stability line. We 
show that the growth of the neutron skin in unstable nuclei does not obey the saturation condition 



because of the polarization effect. The value of the neutron-skin thickness Ar„p = \/lr^ — y (r^) 
is caused by the different radii (skin effect) and only slightly by the different shapes (halo effect) 
of neutron and proton distributions. The relative contribution of both effects depends on the 
competition between the symmetry energy, and the spin-orbit and Coulomb interactions. The 
calculations of the isovector shift of the nuclear radius Ar„p show its primarily linear dependence 
on the asymmetry parameter X = {N — Z)/A. 
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I. INTRODUCTION 



Our knowledge about the properties of neutron excess in heavy nuclei and its relation to 
the neutron-rich nuclear matter and the isotopic symmetry energy is still strongly limited. 
In heavy stable nuclei, the average changes in binding energy E and nuclear radius R with 
nucleon content obey the saturation properties. The volume part i?voi of binding energy 
and the nuclear volume itself are proportional to the particle number A with E^oi = — by A 
and R = r^A^^^, where by > and tq are the constants. Both values of by and tq depend, 
however, on the isotopic asymmetry parameter X = {N — Z)/{N + Z). This is because 
of the difference in saturation bulk density, po ~ i^q ^, of nuclei with different values of X. 
The saturation density po is smaller beyond the beta-stability line for neutron-rich nuclei 
where more neutrons are pushed off to form the "neutron coating". One can expect then 
that the growth of the neutron skin in neutron-rich nuclei violates the saturation property 
R ~ A^^^ for the nuclear radius providing a relative shift of both neutron and proton 
distributions [ij. The main characteristic of the neutron skin is the neutron-skin thickness 



^T"np = \^^J^ — Y i^p)^ where -\/(rf) and y (r^) are the neutron and proton root mean 
square (rms) radii , respectively. The value of Ar„p can be caused by the different radii (skin 
effect) and the different shapes (halo effect) of neutron and proton distributions; see also 
Refs. 2|-l6|. The relative contribution of both effects depends on the competition between 
symmetry energy, spin-orbit and Coulomb interactions 7, ^. 

In the present paper we study a deviation of nucleon distribution from the saturation 
behavior in neutron-rich nuclei. We consider the influence of the spin-orbit and Coulomb 



forces on the neutron, \/(r^, and proton, y{rp'), rms radii as well as the relation of 
the shift Ar„p to the surface symmetry energy. We study also the related problems of 
the nucleon redistribution within the surface region (nuclear periphery), in particular, the 
neutron coating and the neutron excess for the nuclei far away from the /3-stability line. 

We combine the extended Thomas-Fermi (ETF) approximation which takes into consid- 
eration the corrections up to the order of h"^ and the direct variational method assuming that 
the proton and neutron distributions are sharp enough, i.e., that the corresponding densities 
Pp(r) and p„(r) fall from their bulk values to zero in a thin region around the surface. In 
our consideration, the thin-skinned densities Pp(r) and Pn(i") are generated by the profile 
functions which are eliminated by the requirement that the energy of the nucleus should 



be stationary with respect to variations of these profiles. Note that the use of the direct 
variational method and the trial profile function for the particle density allows us to derive 
the equation of state (the dependence of the pressure on the bulk density) in the case of the 
finite diffuse layer of the particle distribution in finite nuclei. 

This paper is organized as follows. In Sec. II we formulate the direct variational principle 
for the density profile function within the extended Thomas-Fermi approximation. Using 
the leptodermous assumption, we obtain A^^^ expansion for the symmetry energy and the 
related values. The results of numerical calculations are presented in Sec. III. We conclude 
and summarize in Sec. IV. 



II. DIRECT VARIATIONAL APPROACH 

We will use the extended Thomas-Fermi approximation which is one of the practical 
realizations of the general Hohenberg-Kohn theorem [9I] on the unique functional relation 
between the ground-state energy and the local density of particles for any fermion system. 
The key point of the ETF is that the total kinetic energy of the many-body fermion system 



is given by the semiclassical expression 10 



12| as follows: 



El,in{Pn,Pp} = -Ekin{Pg, Vpg} = j Ckin [Pn (r) , Pp(r)] , 



(1) 



where ekm [Pn,Pp] = ekm,n[Pn] + ekin,p[Pp], and 
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Here Pq is the nucleon density with g = n for neutron and q = p for proton. 



(2) 



'he semiclassical 



consideration gives the value of parameter (3 in Eq. (|2]) (3 = 1/36 |10|, We point out 
that in the asymptotic limit r — )■ 00, the semiclassical particle density Pq with /3 = 1/36 
goes significantly faster to zero than the one from the quantum-mechanical calculation. 
An asymptotic solution for the particle density pETF,g(i") within the semiclassical extended 
Thomas- Fermi approximation in the limit r — )■ 00 has the following form 11| : 



Petf,, 



exp 



2m Xq 



(3) 
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where Ag is the chemical potential which is negative for a bound Fermi system. A quantum- 
mechanical wave function whose bound s orbital has energy e, gives the partial contribution 
to the asymptotic particle density, which is 



(4) 

Thus, the semiclassical particle density pETF.glr) of Eq. ([3]) for /3 = 1/36 goes faster 
to zero than the quantum- mechanical one (jl]). To overcome this defect of the extended 
Thomas-Fermi approximation the value (3 can be considered as an adjustable parameter. 
We will apply both the semiclassical value /3 = 1/36 and the phenomenological one (3 = 1/9 
which is consistent with the quantum- mechanical asymptotic behavior given by Eq. (jll). 

We will follow the concept of the effective nucleon-nucleon interaction using the Skyrme- 
type force. The functional of the total energy of charged nucleus is given by 

Etot{pg, VpJ = E],in{pg, VpJ + ^pot{Pg, VpJ + ^c{Pp}, (5) 

where Epot{Pq,^Pg} is the potential energy of A^A^ interaction 

^pot{Pq, VpJ = y rfr epot[Pn(r),Pp(r)], (6) 

Cpot [Pn (r), Pp(r)] is the density of the potential energy of the nucleon-nucleon interaction and 
Ec{pp} is the Coulomb energy. In our consideration, the potential energy Epot{Pq,^ Pq} 
includes the energy of the spin-orbit interaction also. 

Following the direct variational method, we have to choose the trial function for Pq{r). 
We will assume a power of the Fermi function for Pg(r) as 

-V 

(7) 

where po,g, Rq, ag, and rj are the unknown variational parameters. Considering the asym- 
metric nuclei with X = {N — Z)/A ^ 1, we will introduce the isotopic particle densities, 
namely the total density p+ = p„ + pp and the neutron excess density p- = pn — Pp with 
P- ^ p+. Assuming a small deviation of the isoscalar bulk density po,+ = po,n + Po,p, the 
radii Rg, and the diffuseness parameters ag with respect to the corresponding average values 
of po, R, and a, we introduce the density profile functions p+(r) and p_(r) to be given by 

4 



Ppartirj 



1 

— exp 



2 m 



er 



Po,q 



1 + exp 



Rn 



p+{r) = pof{r)--pi^^ 



Here, 



+ A, 

a 



fir) 



1 + exp 



P-ir) = ^i^W~2^o^;7- 



r-R 



Ak + A, 

a 

(8)" 



(9) 



the values po and pi are related to the bulk density, R is the nuclear radius, a is the diffuseness 
parameter, and A^ = R^ — Rp and A^ the parameters of the neutron skin. The 

profile functions p+{r) and p^{r) have to obey the condition that the number of neutrons 
and protons is conserved. For the ground state of the nucleus, the unknown parameters po, 
Pi, R, a, Ar, Aa, and r] and the total energy E^ot itself can be derived from the variational 
principle 



6{E - XnN - XpZ) = 0, 



(10) 



where the variation with respect to all possible small changes of po, pi, R, a, Aji, Aa, and rj 
is assumed. The Lagrange multipliers and Xp are the chemical potentials of the neutrons 
and the protons, respectively, and both of them are fixed by the condition that the number 
of particles is conserved. 

We point out that the bare nucleon mass m, but not the effective one m*, enters the 
kinetic energy density ekm,g[Pg] in Eq. ([2]) and accompanies the direct variational procedure 
of Eq. ffTOj) . This is due to the fact that the expression (|2]) is directly obtained as a result of 
a Wigner transform to the quantum-mechanical kinetic energy of many-body systems where 
only the bare mass m is available; see, e.g., Ref. lOj. Usually, the effective mass m* appears 
in the self-consistent Euler equations for the particle density pq within the Thomas-Fermi 
approach (or for the single-particle wave functions in the case of Hartree-Fock theory) after 
the implementation of the corresponding variational procedure. The effective mass appears 
there because the part of the self-consistent mean field, which is caused by the non-local 
interparticle interaction, is associated with the single-particle kinetic energy. In contrast, 
in our direct variational method we do not use the Euler equations for the particle density. 
Thereby, the effective mass m* cannot be included in the kinetic energy density ekin,g[Pg] of 
Eq. ([2]) to avoid a twofold account of contributions from the ti— and t2— components of 
Skyrme forces which enter already the potential energy density epot[Pg('")] in Eq. (Q. This 
argument is also applied to the spin-orbit interaction in a finite system where the spin-orbit 



contribution to the energy functional Epot{Pq,^ Pq} of Eq. ([6]) is involved in the direct 
variational procedure as well. 

We will also assume that the leptodermous condition a/R <^ 1 is fulfilled. The total 



energy (equation of state) ([5]) takes then the following form [l3|: 



Etot(po,^)M = eo(po) + bsipo)A-'^^ + [6v,sym(po) + &s,sym(po) A''/'] + Ec{po,X)/A, 

(11) 

where eo(po) is the energy per nucleon of symmetric nuclear matter, bs{po) is the surface 
energy coefficient, &v,sym(Po) is the volume part of symmetry energy coefficient, &s,sym(Po) 
is the surface part of the symmetry energy coefficient, and Ec{po,X) is the total Coulomb 
energy 

EciPo, X) = acipo) (1 - Xf A'/' + OiA^/% a^Po) = (^) ■ (12) 



The equation of state in the form of Eq. ( jTTl) implies that the total energy per particle 
Etot{po,X)/A is minimized with respect to the independent parameters a, Ar, Aa, and rj 
for arbitrary values of po and X. 

The structure of the equation of state (EOS) given by Eq. ( |TT|) is similar to the semiem- 
pirical mass formula which describes the average changes in nuclear binding energy with 
the mass number. However, in contrast to the mass formula, the bulk density po and the 
asymmetry parameter are not necessarily at equilibrium. The symmetry term ~ X"^ includes 
both the volume, &v,sym(Po); and the surface, &s,sym(Po)) contributions. The surface symme- 
try term &s,sym(po) A^^^^X'^ appears in the advanced mass formula by Myers and Swiatecki 
[isl l and it is currently employed in the description of surface properties and isovector 
excitations in finite nuclei; see, e.g., Refs. [igI, 17 |. 



For a given bulk density po, one can derive the beta-stability line X = X*{A,po) by the 
condition 

dEtot{po,X)/A 



= 0. (13) 

A, X=X* 



dX 

Near the beta-stability line, the total energy per particle ffTTj) is written up to the order 

(X - X*)2 as 

Etot(P0,X)/A = Etot(p0,X*)M+ [6v,sym(Po) + &S,sym(Po) A'^/^ - ^^(Po)^'/'] {X - X*f . 

(14) 
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[Note that Eq. (fl^ is written for A = constant.] The energy of the ground state for a given 
value of mass number A is obtained from the additional equilibrium condition 

d 



-E,,,{po,X*)/A 



= 0, (15) 

^,PO=PO,eq 



dpo 

where po,eq is the equilibrium bulk density. 

The parameters Aji and in the profile functions of Eq. (|H]) derive the neutron skin, 
Avnp, and the neutron excess, Ns, in the surface region of the nucleus ("neutron coat- 
ing"). Substituting Eqs. ([H]) and (jO]) into the conservation particle condition and using the 
leptodermous expansion, we obtain for the neutron excess N — Z the following expression 



where 



N - Z ^ Nv + Ns, (16) 
Ny ^ ^R' (^1 + 3Koiv)^ + QkM (I)') pi, (17) 



.R 



Po 
2 ' 



(18) 



+ A, i^Koir]) + 4:Ki{r])^ + 3^2(^7) (- 
where kjIt]) are the generalized Fermi integrals derived in Ref. 13 | 

/•oo 

Kj{r]) = dx x^ [(1 + e")-" - (-1)^' (1 - (1 + e-")-")] . (19) 
Jo 

The first term A'^v' ~ on the right hand side of Eq. ( [T6|) is caused by the redistribution of 
the neutron excess within the nuclear volume while the second one Ns ~ R^ is the neutron 
coating. 

Note that the variational conditions of Eq. ( ITOll leads to an additional dependence of 
the variational parameters Pq, pi, R, a, A^, A^ and rj on the external parameters A and 
X. The values of A^ and A^ depends slightly on the Skyrme force parametrization. In the 
case of the SkM forces we have evaluated the dependence of A/j and A^ on X for A = 120 
numerically and fitted it by the following formula 

Ar{X) ^ 1.34 X + 0.07 fm, A^(X) ^ 0.36 X + 0.53 X^ fm . (20) 

In general, the change of the radius R of the nucleon distribution with the nucleon 
number A is caused by two factors. There is a simple geometrical change of R because of 



R oc y4^/^. An additional change can occur due to the polarization effect (the bulk density 
distortion) with moving away the beta-stability line. In particular, the size of the neutron 
skin is sensitive to the symmetry and Coulomb energies. To see that we expand the total 
energy Etotipo^^)/^ around the saturation density po.™- By keeping only terms quadratic 
in 6po = Po — Po,cq we rewrite equilibrium Eq. (fT4|) as [18] 



E,,t{po,X)/A = Etot(po,eq,^*)M + 



18po,eq 



(Po - Po, 



+ ■ 



A, sym , 
Po,eq 



(X-X*)^(pO-po,eq), (21) 



where Ka is the incompressibility of finite nucleus 

Ka 



2 d'EUPo,X*)/A 

Po,eq 



and PA,sym is the partial pressure related to the symmetry and Coulomb energies 

d 



pIo.^ N,sym(Po) + &S,sym(Po) A - ac{po)A^/^] 

• ^ cpo 



(22) 



(23) 



Pa, sym 

A,pO=PO,cq 

As seen from Eq. ( l2Ti) . a deviation from the beta-stability line (X ^ X*) implies the change 
of the bulk density po. The corresponding change of po is dependent on the incompressibility 
Ka and the partial pressure PA,sym- For an arbitrary fixed value of X, the equilibrium density 
Po,x is derived by the condition 



^E,,,{po,X)/A 

Opo 



= 0. (24) 

A,po=po,x 

Using Eqs. (I2T1) and (!24|) . we obtain the expression for the shift of the bulk density (polar- 
ization effect) in the neutron rich nuclei 

. Pa, sym 



PO,X — P0,eq — 9- 



K 



-{X-X 
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A 



(25) 



In Fig. [U we have plotted the partial pressure -PA,sym(Po) versus the bulk density po 
(partial equation of state) for the nucleus ^^°S'n. The partial contributions to -PA,sym(po) 
from the symmetry volume ~ f^&v,sym(Po)/f^PO) the sjmametry surface ~ dbs,sym{po)/dpo, 
and the Coulomb ~ dac{po)/dpo terms are also plotted in Fig. [H The dashed vertical line 
shows the position po/po,eq = 0.62 of the spinodal instability border where Ka = 0. On the 
left side of this line the nucleus is unstable with respect to the bulk density variations. 

As seen from Fig. [U the equilibrium partial pressure -P4,sym(Po,eq) is positive and thereby 



20| . We point out that in general the sign of the equilibrium 



Po,x < Po,eq', see also Refs. 

partial pressure -P4,sym(Po,eq) depends on the Skyrme force parametrization and this fact can 



be used to fit the Skyrme forces 
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FIG. 1. The partial pressure -PA,sym foi' the nucleus Sn calculated for the SkM parametrization 
of the Skyrme force. The equilibrium bulk density /Oo,e<j — 0.153 fm~^. The dashed vertical line is 
the spinodal instability border. 

III. RADII OF NUCLEON DISTRIBUTIONS AND NEUTRON SKIN 

As above noted, the bulk density po,x is smaller for neutron-rich nuclei; more neutrons 
should be pushed off to enrich the skin providing a polarization effect. The nuclear rms 
radius 

drr^ p^{r)/ j dr p+{r). (26) 



does not necessarily obey then the saturation condition having that a/ (r^) is nonproportional 
to A^/^. As a consequence, the nuclei with significant excess of neutrons exhibit neutron 
coating, i.e., are characterized by larger radii for the neutron than for proton distributions. 
The interest in the neutron coating was recently raised because of expectations that an 
analysis of the neutron coating could permit extrapolating the nuclear properties to neutron 
matter [fQ^. From the point of view of study of the EOS, the coating size could provide 
information on the derivative of the symmetry energy with respect to the particle density 
23| . In general, the neutron coating Ns of Eq. ( IT7|) can indicate the possibility of a 
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FIG. 2. The rms radius of nuclei near the beta-stabihty hne. The thin sohd hne is for the 
beta-stabihty hne, the thick sohd hne is beyond the beta-stabihty hne for three nuclei, and the 
dotted curve is for the neutron drip line. The dashed lines are the rms radii calculated with the 
step distribution ()27p . The calculations have been performed for the SkM parametrization of the 
Skyrme force. 

giant neutron halo which grows with moving away from the beta stability line [l| . 



In Fig. [2] we have plotted (see the thick solid lines) the nuclear rms radii yJr^ obtained 

from Eq. (l26ll for three nuclei. The results of Fig. [2] are only slightly sensitive to a small 
variation of the diffuse layer and we have here assumed that Aa = 0. The thin solid line of 
Fig. [2] represents the rms radius ~ A^/^ along the beta-stability line X = X*{A) which is 
parametrized by X*{A) = 0.17^2/3/(26.5 - 25.6^-^/3 + o.lTA^/^) [24]. The deviation of 
the rms radii (thick solid lines) from the saturation behavior ~ A^^^ (thin solid line) can 
not be related directly to the appearance of the giant neutron halo at the approach to the 
drip line (dotted curve) because we have here assumed = 0. As above noted, there are 
two sources for the change of the radius of nucleon distribution with the nucleon number 
A. The first one is due to a simple geometrical reason and the second one is because of the 
polarization effect; see Eq. f l25|) . To extract a simple geometrical change of the rms radius 
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a/ (r^) we will perform the calculations of a/ (r^) with a step nucleon distribution 



p(r) = poQl?^ - 



(27) 



Then the geometric rms radius calculated with the step function is given by 



gcom V O 



We will normalize the "geometrical" rms radius to the one \l (r^) 



^fiflR*, where 



geom 



i?* is the nuclear radius on the beta-stability line which obeys the saturation behavior 
R* = tqA^/^. Finally we obtain 



The results of calculations by use of Eq. ( !28|) are shown in Fig. [2] with the dashed lines. 
As one can see these results for the nuclei ^^Na, ^^°Sn, and ^osp]^ g^j^g very close to the ones 
on the beta-stability line (thin solid line). The difference between the dashed lines and the 
thick solid ones represents the magnitude of the polarization effect given by Eq. fl25|) . Thus, 



regions of medium and heavy nuclei is caused by the polarization effect which perturbs the 
distribution of the neutron excess. 

To check the origin of the polarization effect, we have plotted in Fig. [3] the relative shift of 
the neutron and proton radii ^r/ R and the diffuseness parameters A^/a versus the neutron 
excess N — Z for the fixed Z = 50. Note that the direct use of the variational procedure 
ffTOj) with the modified profile functions (r) and p_ (r) of Eq. ([8]) is badly converged with 
respect to the variations of the parameter A^. To overcome this difficulty in Fig. [3l we 
have performed the variational calculations by use of the basic trial functions Eq. (I7l). One 
can see from Fig. [3] that the parameter Aa/a = (a„ — ap)/a, i.e., the parameter of the 
shape distribution, is appreciably growing with the growing of the neutron coating whereas 
the skin parameter Ar/R is only slightly sensitive to the increase of N — Z. We point out 
also that both shift parameters Ar and A^ are only slightly sensitive to the semiclassical 
gradient parameter /3 in Eq. ([2]). 

The sensitivity of the rms radii of the nucleon distribution 




(28) 



we can conclude that the deviation of a/ {r"^) from the saturation behavior ~ A^^^ in the 




(29) 
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FIG. 3. Dependence of the relative shift of the neutron and proton radii Aj^/R = {Rn — Rp)/R and 
the diffuseness parameters Aq/o = {an—ap)/a versus the neutron excess N—Z for the isotopes with 
Z = 50. The dotted hne indicates no spin-orbit interaction and the dash-dotted hne no Coulomb 
interaction for the SkM parametrization. The dashed line is for the SLy230b parametrization. The 
vertical dotted line is for the neutron drip line. 

to the structure of the interparticle interaction is shown in Figs. [Hand [5] for Na isotopes. We 
can see from Fig. Hlthat the ETFA results for a/ (r^) agree quite well with the experimental 
data from The sensitivity of the calculation of y^(rfy to the choice of the Skyrme forces 
for two parametrizations, SkM and SLy230b, can also be seen. Such kind of sensitivity can be 
used to fit the Skyrme force parameters. The two additional lines in Fig. Hlshow the influence 
of the spin-orbit and Coulomb interactions on a/ (r^). As was mentioned above, the spin- 
orbit interaction leads to the deeper potential in the surface region and therefore a nuclear 
core attracts external coating neutrons. This effect reduces the value of the polarization 
effect. The Coulomb interaction acts in the opposite direction. The last is because the 
Coulomb interaction increases the mean distance between nucleons and thereby the rms 

12 




FIG. 4. The rms radius of neutron distribution in Na isotopes for the SkM parametrization (soUd 
hne). The dotted hne indicates no spin-orbit interaction and the dash-dotted hne no Coulomb 
interaction for the SkM parametrization. The dashed hne is for the SLy230b parametrization. 



radius. 



The analogous results for the charge rms radius \ {fp) of Na isotopes are shown in Fig, 



O A small increase of rms radius y {r^^ of the proton distribution with an increase of the 
neutron number is caused by the neutron-proton attraction. Note that the experimental 



data for proton rms radius y (r^) in Fig. |5] manifests the non-monotonic behavior which 
is due to the shell effects. The influence of the spin-orbit and Coulomb interactions on the 
charge radius is the same as in the previously observed case for the rms neutron radius. 
Note also that such kind of behavior of J (r^) is correlated with the A dependence of the 



nuclear Coulomb radius Rc] see Ref. 

As mentioned in Sect. 1, the value of isotopic shift of radii Ar„p can be caused by both 
the skin effect and the halo effect in the neutron and proton distributions. To separate these 
effects, we will represent the value of Ar„p as 

AVnp = Arnp,R + Arnp,a, (30) 
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FIG. 5. The same as in Fig. 4 but for proton distribution. 

where Arnp,R and Arnp,a are caused by the different radii (skin effect) and the different 
diffuseness (halo effect) of the neutron and proton distributions, respectively. The corre- 
sponding values are given by 

Ar^p,n^^[l + l[4iv)~2K,{v)] {^f] ^R, (31) 

and 

Ar„p,a ^ yi [ko{v) - 7 [4(7]) - 2k,{7])] A,, (32) 

In Fig. [6] we have plotted the values of Arnp,R and Arnp,a versus the neutron excess for 
Sn isotopes. 

As can be seen in Fig. [6l the relative contribution of the shape (halo) effect, i.e., Ar„p (j, 
to the isotopic shift of radii Ar^p depends strongly on the parameter /3 of the diffuse tail in 
the nucleon density distribution; see Eq. (I3l). For the semiclassical value of /3 = 1/36, the 
halo effect is quite small near the beta stability line and can play an appreciable role close to 
the drip line only. The situation is significantly different in the case of the phenomenological 
value of /3 = 1/9 where the contribution of Ar„pa, and thereby the halo effect, to Ar„p is 
more appreciable. 
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FIG. 6. The partial contribution to the isotopic shift of radii Ar„p from the skin effect, Ar„p^/j/Ar„p 
(soHd hnes), and the halo effect, Ar„p^a/Ar„p (dashed lines), versus the neutron excess in Sn 
isotopes. The calculations have been performed for two values of the gradient parameter /3 = 1/36 
and /3 = 1/9 indicated near the corresponding lines. The vertical dotted line is for the neutron 
drip line. 



The A dependence of the size of the neutron coating Ar^p = a/ (r^) — y (r^) is illustrated 
in Figs. U\ El and [9] for Na, Sn, and Pb isotopes, respectively. The numerical results have 
been obtained from Eq. (!29l) by use of the basic trial functions Eq. (I7l). The experimental 



data have been taken from Refs. 
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29|. 



As can be seen from Figs. U\ El and [9l the Coulomb interaction affects the isovector shift 
of nuclear radii weakly but with growing of A and X this influence slightly increases. The 
last is because the Coulomb interaction increases the distance between protons, i.e., (rp), 
and reduces thereby the isovector shift. The spin-orbit interaction produces the same effect 
as the Coulomb interaction but with stronger magnitude. As was mentioned above, the 
spin-orbit interaction leads to a deeper potential near the surface region and the nuclear 
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FIG. 7. Isovector shift of nuclear rms radius Ar^p = \J (r^) — y (^r^) in Na isotopes for the SkM 
parametrization. Sohd hne 1 was obtained by use of trial functions Eq. ([7|), i.e., Aq 7^ 0, and 
/3 = 1/36. The dotted line is the same but without spin-orbit interaction and the dash-dotted line 
is without Coulomb interaction. The solid line 2 is for A^ = and /3 = 1/36; the solid line 3 is for 
Aa = and /3 = 1/9. 

core attracts the external neutrons decreasing the diffuse layer of the neutron distribution. 
That reduces the isovector shift of nuclear radii because of Eqs. f l30|) and fl32|) . The spin- 
orbit effect on Ar„p increases with X because the increase of X leads to the contribution to 
the density Pn(^) of neutrons with higher angular momentum. 

The value of the halo effect in Ar„p can be estimated from Figs. [71 El and|9]by comparison 
of the solid lines 1 and 2. The curve 2 was obtained neglecting the contribution from the 
isovector diffuseness distortion, i.e., for = 0. We can see that the halo effect, which occurs 
because of Aa 7^ 0, is quite small even for the lighter nucleus Na. As seen in Figs. [71 [HI and 
[9l the halo effect leads to an increase of the isovector shift Ar„p of nuclear radii. The change 
of the slope of curves Ar„p(yl) due to the halo effect can be also obtained analytically from 
Eqs. ( 130|) . ( 13T|) . and ([32]) . Taking into account a typical dependency of the values Aj:j(X) 
and Aa(X) on the asymmetry parameter X, see e.g. Eq. (1201). and the fact that /to!'?) = ~1 
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and Ki{ri) ~ 1.65 for rj = 2, one can see that the slope of curves Ar„p(yl) grows with an 
increase of the diffuse layer shift A„. Comparing the relative shifts of curves 1 and 2 in 
Figs. [71 [HI and [HI one can also conclude that the halo effect is reduced strongly for the 
heavier nuclei. This is because of a general increase of the nuclear stiffness with respect to 
the variation of the diffuse layer with growing A. 

A growth of the parameter /3 which is responsible for the diffuse tail in the nucleon density 
distribution, see Eq. (|3]), increases the slope of the curve Ar„p(yl); see lines 2 and 3 in Figs. 
[71 [HI and M This is caused by the fact that the diffuse layer of the neutrons exceeds the one 
for the protons and the multiplication of the nucleon density Pq{r) by a factor of in Eq. 
(1291) leads to an emphasis of the peripheral region of the particle density which is stronger 
for the neutron peripheral region than for the proton one. Note also that the crossing point 
in Fig. [7] for isotopes of Na happens for Ar„p = where N = Z. 

IV. SUMMARY 

We have applied the direct variational method within the extended Thomas-Fermi ap- 
proximation with effective Skyrme-like forces to the description of the radii of nucleon dis- 
tributions. In our consideration, the thin-skinned nucleon densities Pp(r) and Pn(r) are 
generated by the profile functions which are eliminated by the requirement that the energy 
of the nucleus should be stationary with respect to variations of these profiles. An advan- 
tage of the used direct variational method is the possibility to derive the equation of state 
for finite nuclei: dependence of the binding energy per particle or the pressure on the bulk 
density Pq. We have evaluated the partial pressure P4_sym which includes the contributions 
from the symmetry and Coulomb energies. The pressure PA,sym is positive driving off the 
neutrons in neutron-rich nuclei to the skin. 

Using the leptodermous properties of the profile nucleon densities Pp(r) and p„(r), we 
have established the presence of the neutron coating Ns- The size of the neutron coating is 
growing with moving away from the beta stability line. In Fig. [2]this fact is demonstrated as 
a deviation of the rms radius of the nucleon distribution from the saturation behavior ~ A^^^ 
in the nuclei beyond the beta-stability line. Moreover, the neutron skin develops by diffusing 
the neutron surface against the changeless proton diffuseness and can be responsible for the 
giant neutron halo in neutron-rich nuclei. 
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The average behavior of the nucleon distribution y (r^^ and the size of the neutron skin is 
satisfactorily described within the extended Thomas- Fermi approximation. The sensitivity 
of the calculations of nuclear rms radii a/ (r^) to the choice of the force parametrization 
can be used to a fit of Skyrme forces. We have pointed that the charge radii of proton 
distributions show the shell oscillations with A which are related to the shell effects in the 
Coulomb energy. The charge radii are connected to the isospin shift of neutron-proton 
chemical potentials AA = A„ — Ap for nuclei beyond the beta-stability line by fixed value of 
the mass number A jl8] . It was shown that the isovector shift of the nuclear radius Ar„p is 
primarily linear dependent on the asymmetry parameter X. The Coulomb and spin-orbit 
interactions do not affect significantly the isovector shift of the nuclear radius. 

We have established the infiuence of the polarization effect given by Eq. fl2Sl) on the 
rms radius ■>/ (r^) of the neutron distribution. This effect increases with the asymmetry 
parameter X and can be responsible for the appearance of the giant neutron halo in the 
nuclei close to the drip line. We have also estimated the relative contribution to the value of 
the isotopic shift of radii Ar„p obtained from both the skin effect and the halo effect. The 
halo effect gives usually a minor contribution to the shift Ar„p and it can be comparable 
with the skin effect near the drip line only. 
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